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Energy partition in low energy fission 
M. Mirea 

Horia Hulubei National Institute for Physics and Nuclear Engineering, 077125 Bucharest, Romania 

The intrinsic excitation energy of fission fragments is dynamically evaluated in terms of the time 
dependent pairing equations. These equations are corroborated with two conditions. One of them 
fixes the number of particles and the another separates the pairing active spaces associated to the 
two fragments in the vicinity of the scission configuration. The fission path is obtained in the frame 
of the macroscopic-microscopic model. The single particle level schemes are obtained within the two 
center Woods-Saxon shell model. It is shown that the available intrinsic dissipated energy is not 
shared proportionally to the masses of the two fission fragments. If the heavy fragment possesses 
nucleon numbers close to the magic ones, the accumulated intrinsic excitation energy is lower than 
that of the light fragment. 
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I. INTRODUCTION 

Under the action of a mutual Coulomb repulsion, at 
scission the fission fragments are accelerated in opposite 
directions. These fragments are highly excited, as under- 
lined in many review papers (T], 0] . The maximal kinetic 
energy issued in the process amounts to the Q-value in 
the case of cold fission. The fragments decay on their 
ground state mainly by evaporation of neutrons and by 
radiation emission. It is known that the motion of any 
physical system is governed by conservative forces and by 
frictional ones that give rise to dissipation. Therefore, the 
excitation energy of the fragments must depend on the 
dynamics of the nuclear system in its path to scission. 

The fission process offers a possibility to investigate 
how two nuclei in contact share their excitation energy. 
In the analyzes of experimental data, the authors of Ref. 
[H evidenced an energy sorting mechanism based on sta- 
tistical arguments. Considering a postulated difference 
between the temperatures of the two nascent fragments 
in conjunction with the condition of maximum entropy, 
they emphasized a flow of energy from one fragment to 
another. This flow of energy depends on the available 
states in the fragments in the scission configuration. In 
this context, an explanation for the violation of the con- 
stant temperature hypothesis Q that involves a propor- 
tionality between the intrinsic excitation energy division 
and the mass ration of the fragment is offered. 

Experimental direct indications about the excitation 
energies of the fragments are obtained by measuring their 
evaporated neutrons 0. Despite a similar temperature 
of the neutron velocity distributions, the experiment re- 
vealed that a larger excitation energy characterizes the 
light mass distribution in comparison with the heavy one. 
The shifting in the sawtoothlike behavior of the neutron 
multiplicity as function of the parent excitation were at- 
tributed mainly to the deformation energy, and not to 
the intrinsic heat. The thermal neutron induced fission 

233u gj analyzed in Ref. Q evidenced a small neutron 
multiplicity in the ^4=132 region. By increasing the exci- 
tation energy of the compound nucleus, it was observed 



that in this mass region the kinetic energy decreases by 
about 2 MeV, that is a large value. The interpretation as- 
cribed a similar temperature of fragments as scission and 
it was speculated that this increment in excitation energy 
is due to a modification of the shape sequences during fis- 
sion leading to a deformed heavy fragment. Furthermore, 
the multiplicity obtained for two neutron induced ^^^Np 
fission energies [8| revealed a modification of the heat of 
only the fragments in the heavier mass distribution. It 
is an experimental indication for a sorting mechanism in 
the intrinsic excitation energy. 

Motivated by these aspects, in this work, the intrinsic 
excitation energy of the fragments are evaluated dynam- 
ically in terms of the time dependent pairing equations 
in the cold fission regime. The macroscopic-microscopic 
is employed to obtain the fission path using the minimal 
action principle. The method is briefly described in the 
next section. The basic ingredients for the time depen- 
dent pairing equations are the single-particle diagrams 
that must be computed from the initial state of the fis- 
sion nucleus up the configuration given by two separated 
fragments. The Woods-Saxon two center shell model @ 
used to determine realistic level scheme along the fission 
path is presented in section lllll In the section ITVl the for- 
malism concerning the time dependent pairing equations 
is introduced and its relevance in calculating dissipation 
energy is emphasized. In section |Vl the formalism is ex- 
tended for two separated nuclei. In Section IVTl the results 
concerning the ^^*IJ fission are reported. The last section 
is devoted to conclusions. 



II. THE FISSION TRAJECTORY 

In order to calculate the energy levels diagrams for the 
fissioning system, the first step is the determination of 
a fission path that satisfies the minimal action criteria 
. The sequence of shapes that follow a nucleus when 
it passes from the ground state to the scission point de- 
pend principally on the potential energy surface and the 
inertia. 
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In the macroscopic-microscopic method, the whole sys- 
tem is characterized by some collective coordinates that 
determine approximately the behavior of many other in- 
trinsic variables. The basic ingredient in such an ana- 
lyzis is the shape parametrization that depends on several 
macroscopic degrees of freedom. The generalized coordi- 
nates associated with these degrees of freedom vary in 
time leading to a split of the nuclear system in two sepa- 
rated fragments. The macroscopic deformation energy is 
calculated within the liquid drop model. A microscopic 
potential must be constructed to be consistent with this 
nuclear shape parametrization. A microscopic correction 
is then evaluated using the Strutinsky procedure. 

First of all, in our description, it is required to define a 
nuclear shape parametrization. In the following, an axial 
symmetric nuclear shape is obtained by smoothly join- 
ing two spheroids of semi-axis and bi {i—1,2) with a 
neck surface generated by the rotation of a circle around 
the axis of symmetry. By imposing the condition of vol- 
ume conservation we are left by five independent gener- 
alized coordinates {qi} (i=l,5) that can be associated to 
five degrees of freedom: the elongation R given by the 
distance between the centers of the spheroids; the neck- 
ing parameter C = S/R3 related to the curvature of the 
neck, the eccentricities associated with the deforma- 
tions of the nascent fragments and the mass asymmetry 
parameter rj = 01/02. The notations that describe this 
parametrization can be identified by inspecting the Fig. 
1. Due to the axial symmetry, the surface equation is 
given in cylindrical coordinates for the three regions in- 
volved: 



p{z) 



P3 - S^Hj - [z - z^ Y 
&2\/l - {z- Z2Y/al, 



z < Zci; 

Zcl < Z < Zc2; 

Z > Zc2. 



(1) 



It is known that a nuclear shape is well adapted for 
the fission process is the following conditions are satis- 
fied [ll|: (i) The three most important degrees of free- 
dom, that is, elongation, necking and mass-asymmetry 
are taken into account; (ii) A single sphere and two sepa- 
rated fragments are allowed configurations; (iii) The flat- 
ness of the neck is an independent variable. All these 
conditions are fulfllled by the above parametrization. If 
5"=!, the shapes are necked in the median surface char- 
acterizing scission shapes and for S—-1 the shapes are 
swollen addressing the ground state and the saddle con- 
figurations. 

If we consider that the elongation qi = R is the main 
coordinate, the dependencies of the other generalized co- 
ordinates Qi = fi(R) {i = 2, 5) must be obtained. As 
specified in Ref. [l^l, such trajectories emerge by mini- 
mizing the action functional. 



^2M{q,,dq^/dR)V{q,)dR (2) 



where M{qi,dqi/dR) is the inertia along the trajectory 
and V{qi) is the deformation energy. Ri and Rf stand 




FIG. 1: Nuclear shape parametrization. Two ellipsoids of 
different eccentricities are smoothly joined with a third sur- 
face. Two cases are obtained: (a) the curvature of the neck 
is positive 5*=! and (b) the curvature is negative S—-1. 



for the elongation associated to the ground state and the 
exit from the barrier, respectively. In our calculation 
the reference of the deformation energy is always taken 
as the energy in the ground state. So the next condi- 
tion is fulfilled V{Ri) = V{Rf) = 0. As it can be seen 
in formula ([2]), as the fissioning nucleus passes from its 
ground state to the scission configuration, the sequences 
of shapes depends mainly on the deformation energy and 
the inertia. The deformation energy is obtained in the 
frame of the macroscopic-microscopic model |12l] while 
the inertia is computed within the cranking approxima- 
tion fl3, [13, • The deformation energy was obtained 
by summing the liquid drop energy El dm with the shell 
and the pairing corrections 5E. 



V = E 



LDM 



SE 



(3) 



The macroscopic energy Eljjj^i is obtained in the frame- 
work of the Yukawa - plus - exponential model (l6l | ex- 
tended for binary systems with different charge densities 
as detailed in Ref. 1171: 



Eldm = En + Ec + Ev 



(4) 
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where 



En 



0-2 



'0" JviJv2 

is the nuclear term, 



(5) 



The total inertia is the sum of the contributions that cor- 
respond to the proton and to the neutron level schemes. 
Usually, the matrix elements of the derivatives of the 
Hamiltonian in Rel. ^ are replaced by the matrix ele- 
ments of the derivatives of the mean field potential alone. 



E, 



c 



ri2 



(6) 



is the Coulomb energy, and Ey is the volume energy. 
In the previous definitions pe are charge densities and 
ri2 = 1^1 — r2|- The numerical values of the parameters 
02, To, and a are taken from Ref. p^ . 

The shell effects 6E are obtained as a sum between the 
shell and the pairing microscopic corrections. In this con- 
text, the Strutinsky procedure [l^ was used. These cor- 
rections represent the varying parts of the total binding 
energy caused by the shell structure. The single parti- 
cle level diagrams are computed within the Woods-Saxon 
superasymmetric two-center shell model. 

The effective mass is com put ed within the cranking 
adiabatic approximation [Tol. Il3l. [T^. In a multidimen- 
sional deformation space, where the nuclear shape is de- 
scribed by the set of n independent generalized coordi- 
nates Qi , the inertia tensor Mij is defined by the equation 
of the kinetic energy T: 



T : 
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In the adiabatic description of the collective behavior of 
a nucleus, the nucleons are assumed to move in a aver- 
age deformed potential. Using a Hamiltonian H{qi, ...qn) 
that includes pairing interactions, introducing the collec- 
tive parameters qi by means of the Lagrange multipliers, 
it is possible to obtain the response of the nuclear system 
for slow changes of the shape within the cranking model 
formula 



My ((71, ...,q, 



\0y 



(8) 



where \v > and |/i > are single particle wave functions. 
El,, Ui, and Vi, are the quasiparticle energy, the vacancy 
and occupation amplitudes of the state h', respectively, in 
the BCS approximation, and Pij is a correction that de- 
pends on the variation of the pairing gap and the Fermi 
energy as function of the deformation coordinates. Re- 
cently, the formula ([S]) was generalized by taking into 
account the intrinsic excitation produced during the fis- 
sion process itself ^15]. The inertia M along a trajectory 
in the configuration space spanned by the generalized co- 
ordinates qi (i=l,5) can be obtained within the formula 



M 



5 5 

i=i j=i 



dqi dqj 
dR dR 



(9) 



III. SINGLE PARTICLE ENERGIES 

A microscopic potential must be constructed to be con- 
sistent within our nuclear shape parametrization. The 
simplest way it to use a semi-phenomenological Woods- 
Saxon potential. In order to take into account nuclear 
deformations going over to separate shapes and obtain 
two separated fragments, a two-center shell model with a 
Woods-Saxon potential was developed recently Q. Other 
recipes that allows to treat strongly deformed nuclei are 
presented in Ref. [l^ [2^. The mean field potential is 
defined in the frame of the Woods- Saxon model: 



V, 



I + cxp 



A(p,2) 



(10) 



where A(p, z) represents the distance between a point 
(/9, z) and the nuclear surface. This distance is measured 
only along the normal direction on the surface and it 
is negative if the point (p, z) is located in the interior 
of the nucleus. 14 is the depth of the potential while 
a is the diffuseness parameter. In our work, the depth 
is Vc = Voc[l ± K(iVo - Zo)/Nq + Zq)] with plus sign 
for protons and minus sign for neutrons, Vbc= 51 MeV, 
a=0.67 fm, k=0.67. Here ^o, -^o and represent the 
mass number, the neutron number and the charge num- 
ber of the parent, respectively. This parametrization, 
referred as the Blomqvist-Walhlborn one in Ref. [2l|, is 
adopted because it provides the same radius constant ro 
for the mean field and the pairing field. That ensures a 
consistency of the shapes of the two fields at hyperde- 
formations, i.e., two tangent ellipsoids. The Hamiltonian 
is obtained by adding the spin-orbit and the Coulomb 
terms to the Woods-Saxon potential. The eigenvalues 
are obtained by diagonalization of the Hamiltonian in the 
semi-symmetric harmonic two center basis [23 - l23 |. In 



this work, the major quantum number used is Nmax=^'^- 
The two center Woods-Saxon model will be used to com- 
pute shell and pairing corrections together with inertia 
in this work. The two center shell model represents a 
valuable instrument to investigate the role of individual 
orbitals for the treatment of a wide variety of nuclear pro- 
cesses like cold fission (2^ , formation of superheavy ele- 
ments [26[ or superasymmetric disintegration processes, 
pertaining to cluster- and alpha-decays [27l - [29| . 



IV. TIME DEPENDENT PAIRING EQUATIONS 

In the actual formalism, the starting point is a many- 
body Hamiltonian with pairing residual interactions. 
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This Hamihonian depends on some time-dependent col- 
lective parameters q{t) = {qi{t)} (i = 1, ...n), such as the 
inter-nuclear distances between atoms or nuclei: 



k>0 



kA>0 



(11) 

Here, ek are single-particle energies of the molecular po- 
tential, and ak denote operators for creating and de- 
stroying a particle in the state k, respectively. The state 
characterized by a bar signifies the time-reversed part- 
ner of a pair. The pairing correlation arise from the 
short range interaction between fermions moving in time- 
reversed orbits. The essential feature of the pairing cor- 
rection can be described in terms of a constant pairing 
interaction G acting between a given number of particles. 
In this paper, the sum over pairs generally runs over the 
index k. Because the pairing equations diverge for an 
infinite number of levels, a limited number of levels are 
used in the calculation, that is N levels above and below 
the Fermi energy Ep. 

In order to obtain the equations of motion, we shall 
start from the variational principle taking the following 
energy functional 



d 

C —< (p \ H — ifi-'o^ I ^ > 



(12) 



by assuming the many-body state as time dependent BCS 
seniority-zero wave function 



'Pit) > = '[[{uk{t)+Vk{t)a 



k'^k. 



(13) 



To minimize this functional, the expression (|12l) is de- 
rived with respect the independent variables Vk , together 
with their complex conjugates, and the resulting equa- 
tions are set to zero. Eventually, the next coupled- 
channel equations are obtained : 



ihpk 
ihkk 



KkA* 

i2pk 



1)A - 2Kkek 



(14) 



where pfc =| Vk |^ are occupation probabilities, Kk = u*j,Vk 
are pairing moment components, and IS. — G Kk is the 
pairing gap. Uk and Vk are the complex BCS occupation 
and vacancy amplitudes. The variations of single-particle 
densities pk can be evaluated for different values of the 
generalized velocities by solving the previous system of 
coupled equations as done in Ref [30,]. Eqs. (|T4l) are also 
generically known as the time dependent Hartree-Fock- 
Bogoliubov equations [sil,!!^]. As mentioned in Ref. 32], 
a connection with the Landau-Zener effect is included in 
these equations. Levels undergo Landau-Zener transi- 
tions on virtual levels with coupling strengths given by 
the magnitude of the gap A. Recently, these equations 
were generalized to take into account the Landau-Zener 
effect in seniority one systems [1, and the pair break- 
ing mechanism [34|. 

These time-dependent pairing equations can offer a 
measure of the dissipated energy. The difference between 



the total energy value E obtained within the TDHFB 
equations and Eq given by the static BCS-equations rep- 
resents an approximate measure for the dissipation Ejj: 

Ed=E- Eo. (15) 

E is expressed simply in terms of pk and Hk 



E = 2j2ekPk-G\J2^k \' -gJ2pI 



(16) 



Eq corresponds to pi and k° associated to the lower en- 
ergy state, that is, obtained from BCS equations. 

The sum of single particle densities derivatives of Eqs. 
(fll is 



ifiJ2kPk = Efe KA* - K*A] 
= ^[\A\^-\A\^]=0 



(17) 



This result shows that the sum of the single particle occu- 
pation probabilities is a constant quantity as the system 
evolves in time according to Eqs. (|14p . that is, the av- 
erage number of particles in the pairing active space is a 
constant of the motion. 



V. NUMBER OF PARTICLES 

After scission, the levels from the pairing active space 
will be shared between the two fission fragments. The 
levels of the core are sorted accordingly. The sum of oc- 
cupation probabilities of the levels located in each frag- 
ment must be equal to the number of nucleons. This is 
a problem that can be solved by appealing to two prop- 
erties of the nuclear system. 

First of all, the sorting is produced in a continuous 
manner: the wave function associated to one single- 
particle energy level is transferred gradually to one of 
the two potential wells obtained asymptotically, after the 
scission. For example, the Woods-Saxon wave function 
of the lowest single particle energy is displayed in Fig. 
[5] for different values of the distance between the cen- 
ters of the fragments, that is, for different shapes of the 
potential. When the two fragments are completely sep- 
arated, the wave function is located in one of the two 
wells. Making use of this property, it is possible to iden- 
tify the final localization of the single particle level even 
when the system behave as one single nucleus if the final 
mass-asymmetry is known. For this purpose we calcu- 
late two quantities Qifc =< ^k{z) \ Q{~z) \ (pk{z) > and 
Q2k =< fkiz) I Qiz) I ipk{z) > where 6 is the Heaviside 
function and tpk is the Woods-Saxon single particle wave 
function of the level k along the axis z. If Qifc > Q2k, 
the wave function is located in the well 1. It is worth 
to note that this precedure fails only in the avoided level 
crossing regions. 

Secondly, the matrix element of the pairing interac- 
tion G is in principle dependent of the overlap of the 
wave functions of the pairs [s^ [s^. As long as a sin- 
gle nucleus is involved, the monopole approximation is 
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FIG. 2: Right panel. For the left axis, the lowest energy Woods-Saxon wave function for the two center model (^0(2) as function 
of z for different values of the distance between the centers is displayed with a thin line. The right axis corresponds to a plot 
with a thick curve of a section in the middle of the neutron Woods-Saxon potential Vb(2) at the same values of the distance 
between centers. In the left panel, a representation of the Woods-Saxon potential Vb in the cylindrical p and z coordinates is 
made. The distances between centers R are 0, 6, 12, 18 and 21 fm from top to bottom. The configurations displayed correspond 
to the minimal ^'^''U fission trajectory. 



considered to perform well [37| and all the values of the 
matrix elements are considered to be equal in the active 
level space. However, after the scission, the matrix ele- 
ments of the pairing interaction between wave function 
belonging to different fragments must be zero. If the 
pairing matrix elements between pairs located in differ- 
ent fragments at scission is zero (Gi2=0), then the energy 
given by Rel. becomes 



2 ek.Pk, - Gi I J2k^ P -Gi E/ci Pfei 
+ Efc, efc2Pfc2 - G2 I Efc, '«fc2 I' -G2 Ek, pI, 

= i?i + E2 

where the sum over ki and k2 means that the levels be- 
long to fragments 1 and 2, respectively. We used an ar- 
row in the previous relation to indicate that the constant 
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value G of the parent nucleus can be also transformed 
in the two values Gi and G2 associated to the two frag- 
ments. On other words, the monopole approximation is 
considered valid separately in each fragment by consid- 
ering constant values of Gi and G2. The relation ([18]) 
shows that the total energy E is decomposed in two frag- 
ment total energies Ei and E2 given by relations of the 
type (|16p if the values of the matrix element of the pairing 
interaction addressing different fragments Gi2=0. 

Within these properties, a simple recipe to fix the av- 
erage number of particles at scission can be elaborated. 
Taking as example the proton level scheme, the condi- 
tions that the sum of occupation probabilities of levels in 
the two wells must be equal with the number of nucleons 
of the fragments can be written as: 

Zp, + Zc, = Zi (19) 

Zp2 + = Z2 

where Zi (i—l,2) are the number of protons in the two 
fragments, Z^ and Zp. stand for the number of protons 
in the core and the number of protons in the pairing ac- 
tive space, respectively. For an initial number of pairs N 
considered in the pairing active space, the values of Z^ 
and Zp-^ are simply obtained by counting the levels given 
by the two center model. The occupations probabilities 
Pki must be obtained from the time dependent pairing 
equations HH). Exploiting the two previous properties, 
the problem to fix the final average number of particle 
within equations is now trivial. After the passage 
of the external saddle point, in the descent to scission, 
we insert the condition in the functional (IT^ and 
we continue to solve the equations of motion. When the 
good numbers of particles are obtained, we impose the 
condition that G12 is zero between the wave functions be- 
longing to separated fragments. Both conditions (IT9|) is 
introduced after the passage of the external saddle point. 
In the operator notation the condition (|T5)) becomes 

Zp2 Zp-^ Zp-^ Zp^ 

= EfeiKi< +«fci4i) (20) 

Zp2 = Y.k2i^k2ak2 + "fcs^t) 

This condition is introduced in the energy functional (jl2p 
d 

£ =< LP \ H - ih— - \{Zp.^Zp^ - Zp^Zp^) I LP > (21) 

using a Lagrange multiplier A. Imposing also the condi- 
tion that the interaction matrix element G between pairs 
of the same fragment is not the same than those of dif- 
ferent fragments the new time dependent equations read, 
eventually: 

ihpk^ = Kk^Al - KfcjAi 

ihpk2 = Kk2^*2 - 'tfc^Aa ^22) 
iMfci = (2pfcj - l)Ai - 2Kfci(efe^ - AZpJ 
ihkk2 = {2pk2 - 1)^2 - 2Kk2{ek2 + A^pJ 



where Ai = Gi J^ki '^fci + ^12 J2k2 '^^2 and A2 = 
G12 J2ki + ^2 J2k2 '^^2- When Gi2=0, it can be eas- 
ily verified the average number of particles in the two 
fragments are conserved according to conditions of the 
type P7|) applied separately on the two working spaces. 
That means, Zi and Z2 behave as constants. The previ- 
ous recipe represents the simplest dynamical method to 
project the average number of particles onto two separate 
nuclei. 



VI. RESULTS 

To obtain the fission trajectory, the action integral ([5]) 
must be minimized in our five-dimensional space. The 
first turning point Ri is obtained by determining the 
ground state configuration while the second one Rf lies 
on the equipotential surface that characterize the exit 
from the outer barrier. That means, i?/ is defined by the 
multidimensional function V{R, G, ei, €2, 'rj)—0. Different 
methods are currently envisaged to obtain the heights 
of the barriers. In static calculations [s^ the immersion 
procedure is extensively used while for dynamical paths 
j39M4ll | the Ritz method is applicable. To minimize the 
action integral we used a numerical method initiated in 
Ref. [12] and used for fissionprocesses in a large range of 
mass asymmetries [s^, I43l - l45j The dependencies C{R), 
ei(i?), e2(-R), and r]{R) were considered as spline func- 
tions of n variables Ck, eifc, ^2k and rjk {k — l,n) that 
are associated to the fixed mesh points Rk located in the 
interval (i?i, Rj). The integral action is transformed in a 
numerical function that depends on the An variables and 
it is minimized numerically. 

Determination of potential energies V and of effective 
masses Mij are very expensive in computing time. For 
the numerical minimization procedure, a large number of 
iterations is required and it is not possible to calculate 
the values of V and My for each iteration. An interpola- 
tion of calculated values of the energy and of the effective 
masses will be used. Therefore, to make the problem 
tractable, first of all, a grid of deformation values was 
fixed in the five-dimensional configuration space: 25 val- 
ues of R between fm and 24 fm, 7 values of the eccentric- 
ities ti between and 0.6,7 values of the ratio rj in the in- 
terval 1 and 1.6, and 23 values of G between -0.115 fm~^ 
and 0.105 fm~^. The pertinent region of deformations for 
the possible trajectories between the ground state and 
the exit point from the fission barrier was spanned. In 
each point of the configuration space, the deformation 
energy and the tensor of inertia was computed. During 
the minimization process, interpolated values of the de- 
formation energy and the inertia were used. Different 
initial values of the generalized coordinates were used as 
input parameters for the minimization. Therefore differ- 
ent local minimum were obtained. The best value was 
selected. Moreover, additional calculation of the action 
integral were performed by slowing varying the general- 
ized coordinates around the best trajectory previously 
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FIG. 3: ^■^''U fission barrier V for a final partition 
'^"^Zr+^^^Te determined along the minimal action trajectory. 
Some particular shapes related to the ground state, the ex- 
tremes of the barrier, the exit point and the scission point 
are inserted in the plot. The distances for the elongation R 
that characterizes the shapes are 4.17, 7.7, 10.43, 12.64, 15.53, 
17.53 and 20.2 fm. 



obtained from the numerical procedure. Among all re- 
sults, the best final trajectory for the least action was 
retained. Such a procedure was used in determining a 
theoretical systematic of fission barrier heights in Ref. 
[46| . In the present work, the trajectory in the configu- 
ration space was modified after the saddle of the outer 
barrier. The generalized coordinates were changed to ob- 
tain at scission a required final configuration as explained 
below. 

The minimization leads to a definite path in the sub- 
barrier region. Its extrapolation to the scission point 
is not unique [i^. We will choose the scission config- 
uration be searching the best candidates that have de- 
formations close to those obtained in the exit point of 
the barrier. We found that the ecceirtricities and mass- 
asymmetry parameters in the exit point of the outer bar- 
rier are consistent with the formation of a pair given by 
the isotopes ^"^Zr and ^^^Te. The ground state configu- 
rations of the fragments were calculated for this purpose. 
By introducing the ground states of the fragments in the 
scission configuration, the excitation due to the deforma- 
tions energy is minimized and we are left only with the 
dissipation. As evidenced from evaluations |48|, the par- 
tition ^°^Zr-|-^'^^Te is one of the larger enought probable 
partitions. Moreover, this final configuration is also in- 
teresting due to the fact that the numbers of nucleons in 
the heavy fragments are close to the magic ones. This 
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FIG. 4: Neutron level diagram for the ^^^U fission with re- 
spect the elongation R. 
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FIG. 5: Same as Fig. |4]for the proton diagram. 



configuration allows also a simple comparison with the 
hypothesis made in Ref. Q . The fission barrier together 
with some particular shapes are plotted in Fig. |3l 

The neutron and proton single particle diagrams are 
calculated along the minimal action trajectory, from the 
ground state of the parent up to the formation of two 
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FIG. 6: Dissipation energy E* for protons (upper panel) 
and neutrons (lower panel) as function of the elongation 7?. 
The calculations are made within Rel. (|15p and Eqs. (|14|l 
without imposing any condition for fixing the number of of 
particles. The dashed line, the dot dashed line and the full 
line correspond to a inter-nuclear velocity dR/dt of 10*, 10^ 
and lO" m/s, respectively. 



separated fragments. These level schemes are plotted in 
Figs, m and [5] In Fig. |31 at i?=0, the parent nucleus 
is spherical. For small deformations the system behavior 
is similar to a Nilsson diagram for prolate deformations. 
At the right of the figure the orbitals of the parent are 
labeled with their spectroscopic notations. The levels be- 
longing to the heavy fragment can be identified long be- 
fore the scission configuration and are plotted with thick 
lines. The heavy fragment is spherical at scission, so its 
levels are bounced in shells. Making use of the fact that 
the heavy fragment becomes spherical after scission, we 
labeled its shells with spectroscopic notations. The lev- 
els of the light fragment are not degenerate due to its 
deformation. In the proton diagram of Fig. [SJ it can be 
observed a smooth decrease of the single particle energies 
after the scission due to the Coulomb mutual polariza- 
tion. The energy slope for the light fragment is larger 
than that for the heavy one. 

It is known that if the projection on a given parti- 
cle number changes the pairing coupling constant G [50| . 
However, for our exploratory investigation, the same G 
will be kept for the parent nucleus and for the fragments. 
After scission, we calculate the matrix elements for each 
fragment separately. The values of G, Gx^G2 are ob- 
tained from the renormalization procedure [lO[ by taking 
into account the number of levels in the pairing active 
space. 



To solve the Eqs. (|T4)) we need the velocity of the gen- 
eralized coordinate. Our calculations will pertain to the 
cold fission regime and are characterized by small val- 
ues of the collective kinetic energy and of the excitations 
energies. These properties can be obtained by selecting 
an appropriate value for the collective velocity. There- 
fore, different constant values of the inter-nuclear veloc- 
ity R ranging from 10'' to 10^ m/s were tested. These 
values can be translated in a time to penetrate the bar- 
rier ranging m the interval [10"^^, lO'^O] s. The time 
for the descent between saddle to scission reaches about 
4x10^^^ s for a velocity of 10^ ni/s, a value considered 
a typical time for scission [2| . Within the semi-adiabatic 
cranking model [lH|, the inertia in the ground state and 
in the asymptotic region of two separated fragments are 
1.26 and 1.208 ^i^/MeV^/fm^, respectively. Within this 
estimation of the mass, the velocities can also be trans- 
lated in a macroscopic kinetic energy that amounts up 
to 0.3 MeV, of the order of magnitude of the zero point 
vibration energy. The system (jl4p is solved numerically 
for the selected velocities. The initial conditions for pk 
and Kk were given by the BCS solutions for the parent 
ground state, a configuration characterized by R «4 fm 
in Fig. [3] In Fig. |6l the dependencies of the dissipated 
energy E* for the different internuclear velocities inves- 
tigated are displayed as function of the distance between 
the centers of fragments R. No conditions for fixing the 
number of particles are used. The dissipated energy for 
the proton level scheme is lower that that of the neu- 
tron one. The dissipated energy for the lower collective 
velocity is negligible, while for the higher one the total 
dissipation reaches about 10.5 MeV. This value is larger 
than in previous estimations [TEj because of the imposed 
modification of the fission trajectory that gives a par- 
ticular configuration at scission. In connection with the 
shape of the barrier displayed in Fig. [3l it can be deduced 
that the larger part of the excitation is formed during the 
penetration of the second barrier and, in the descend be- 
tween saddle to fission. This result is not in line with the 
observation that the excitation energy does not increase 
after the passage of the saddle point 49] . 

It will be interesting to compare the dissipated energy 
and the number of nucleons attributed to the two nuclei 
in the case of the treatment without constraints and the 
treatment with projection of the average number of par- 
ticle. In Fig. El^a) the unconstrained dissipation energy 
for the neutron workspace is plotted with a thick curve 
as function of R at dR/dt— 10^ m/s. The dissipation 
obtained within the scheme proposed for projecting the 
average number of particles is displayed with a thin line. 
Up to the the external saddle, the Eqs. ([14]) are solved . 
After R «16 fm, we imposed the condition (|20p . It can 
be observed that the time dependent pairing equations 
with constraint give a larger dissipation. In Fig. ^h) the 
sums Np^ = J2k, Pki and Np^ = J2k2 Pk^ ^^"^ t^e uncon- 
strained equations are plotted with a thick line, fci and 
k2 are levels in the active pairing space attributed to the 
fragment 1 and the fragment 2, respectively. The total 



number of pairs is conserved Ni + N2 = N = 30. After 
imposing the condition (PO)) . the number of pairs located 
on the two level schemes of the two fragment reach the 
correct values, that is Ni—15 and A^2=15, that define the 
partition analyzed. 

It was remarked in Rcf. 31] that the maxima of the 
dissipated energy arise because of the character of the 
ground state solution is changing rapidly that the sys- 
tem cannot adjust itself. So, the system appears excited 
not because of any changes in the occupation amplitudes 
but because the changes in the ground state to which the 
dynamic state is studied. In this respect, in Fig. [3 we 
analyzed also the behavior of the system by replacing the 
ground state. Instead of the ground state of the parent, 
we used as Eq the sum between the lower energy states 
calculated for the two level schemes of the two nascent 
fragments. The result is plotted with a dashed line. At 
scission we obtain approximately 8 MeV excitation en- 
ergy for both fragments. 

In Fig. [5J the dissipation obtained for the two frag- 
ments are displayed separately. It is found that the dissi- 
pated energy for the heavy fragment is much lower than 
that for the light one. This result is consistent with the 
experimental finding addressing the neutron multiplici- 
ties. It is also revealed in Fig. |S]that in the cold regime 
the increase in excitation energy of the heavy fragment is 
smaller than in the light one, for a modification of the ex- 
citation energy of the parent nucleus. On another hand, 
the experimental findings 'sj shows that the excitation 
energy of the compound nucleus goes almost completely 
into excitation energy of the heavy fragments. That 
means, a large part of this energy is transferred in de- 
formations of the heavy nucleus, not taken into account 
here. A similar behavior is predicted from statistical con- 
siderations in Ref. Q for partitions involving nuclei close 
to magic numbers. 

This model is based on the same philosophical grounds 
as that found in Refs. [13, [5l|. The main differences 
are given by the time dependent equations used and by 
the fact that the dissipation is not evaluated for each 
fragment separately. They used the time dependent 
Schrodinger equations and introduced quasiparticle ex- 
citations through the cranking approximation. 



the dissipation energy in fission fragments is intimately 
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FIG. 7: (a) The dissipated energy without constraints is 
plotted with a full thick line as function of the elongation 
R. The dissipated energy after imposing the condition for 
projecting the number of particles is plotted with a thin line. 
The dissipated energy obtained by replacing the ground state 
of the parent nucleus within the ground states of the two 
nascent fragments, (b) The number of pairs Np^ and Np^ 
located on the levels that belongs to the first and second well 
are plotted with thick lines. The number of particles after 
imposing the conditions for projecting the number of particles 
are plotted with thin lines. Np^ corresponds to the heavy 
fragment and it is always larger than Np^ . 



VII. CONCLUSIONS 

Dynamical estimations of the excitation energies in 
cold fission are evaluated within time dependent pair- 
ing equations. Using conditions that fix the number of 
particles in each fragment it was possible to obtain for 
the first time the excitation energy of each nucleus issued 
in the process. A recent hypothesis that claims the 
excitation energy is not equilibrated between fragments 
was confirmed in the cold fission regime. If the heavy 
fragment is close to magic number, its excitation energy 
is smaller than that of the light one. It was found that 



related on the distribution of pairing occupation proba- 
bilities of the levels at scission. These probabilities can 
be obtained by solving an appropriate set of equations of 
motion. 
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FIG. 8: Excitations energies for "'^Zr (dashed line) and ^^^Te 
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sum between the excitations energies obtained for protons and 
neutrons are presented. 



I.S. Grant, Rep. Prog. Phys. 39, 955 (1976). 
R.W. Hasse, Rep. Prog. Phys. 41, 1027 (1978). 
K.-H. Schmidt, and B. Jurado, Phys. Rev. Lett. 104, 
212501 (2010). 

H.A. Bethe, Phys. Rev. 50, 332 (1936). 

H.R. Bowman, S.G. Thomson, J. CD. Milton, and W.J. 

Swiatecki, Phys. Rev. 126, 2120 (1962). 

V.F. Apalin, Yu.N. Gritsyuk, I.E. Kutinov, V.I. Lebedev 

and L.A. Mikaelian, Nucl. Phys. 71, 553 (1965). 

S.C. Burnett, R.L. Ferguson, F. Plasil, and H.W. 

Schmitt, Phys. Rev. C 8, 2034 (1971). 

A. A. Naqvi, F. Kappeler, F. Dickmann, and R. Miiller, 

Phys. Rev. C 34, 218 (1986). 

M. Mirea, Phys. Rev. C 78, 044618 (2008). 

M. Brack, J. Damgaard, A.S. Jensen, H.C. Pauh, V.M. 

Strutinsky and C.Y. Wong, Rev. Mod. Phys. 44, 320 

(1972). 

U. Brosa, S. Grossman, and A. MuUer, Phys. Rep. 197, 

167 (1990). 

J.R. Nix, Arm. Rev. Nucl. Sci. 22, 65 (1972). 

M. Mirea, R.C. bobulescu, and M. Petre, Rom. Rep. 

Phys. 61 , 646 (2009). 

M. Petre, R.C. Bobulescu, C. Petre, and M. Mirea, Rom. 
J. Phys. 55, 741 (2010). 

M. Mirea, and R.C. Bobulescu, J. Phys. G 37, 055106 
(2010). 

K.T.R. Davies, and J.R. Nix, Phys. Rev. C 14, 1977 



[17 



[18; 

[19 

[2o; 

[21 

[22' 
[23^ 
[24' 
[25| 

[26 

[27 
[28' 
[29' 

[m 

[31 
[32' 
[33' 

[34; 



(1976). 

M. Mirea, O. Bajeat, F. Clapier, F. Ibrahim, A.C. 
Mueller, N. Pauwels and J. Proust, Eur. Phys. J. A 11, 
59 (2001). 

P. MoUer and J.R. Nix, At. Data Nucl. Data Tables 26, 
165 (1981). 

V.V. Samarin, Phys. At. Nucl. 73, 1416 (2010). 

A. Diaz-Torres and W. Scheid, Nucl. Phys. A 757, 373 

(2005). 

S. Cwiok, J. Dudek, W. Nazarewicz, J. Skalski, and T. 

Werner, Comput. Phys. Commun. 46, 379 (1987). 

J. Maruhn, and W. Greiner, Z. Phys. 251, 431 (1972). 

M. Mirea, Phys. Rev. C 54, 302 (1996). 

M. Mirea, Nucl. Phys. A780, 13 (2006). 

M. Mirea, D.S. Delion and A. Sandulescu, Phys. Rev. C 

81, 044317(2010). 

M. Mirea, D.S. Delion and A. Sandulescu, EPL 85, 12001 

(2009). 

M. Mirea, Phys. Rev. C 57, 2484 (1998). 

M. Mirea, Eur. Phys. J. A 4, 335 (1999). 

M. Mirea, Phys. Rev. C 63, 034603 (2001). 

M. Mirea, L. Tassan-Got, C. Stephan and CO. Bacri, 

Nucl. Phys. A 735, 21 (2004). 

S.E. Koonin and J.R. Nix, Phys. Rev. C 13, 209 (1976). 
J. Blocki and H. Flocard, Nucl. Phys. A 273, 45 (1976). 
M. Mirea, Mod. Phys. Lett. A 18, 1809 (2003). 
M. Mirea, Phys. Lett. B 680, 316 (2009). 



11 



[35] D.S. Dclion, M. Baldo, and U. Lombardo, Nucl. Phys. A 

593, 151 (1995). 
[36] R.C. Kennedy, Phys. Rev. 144, 804 (1966). 
[37] J. Dobaczewski, W. Nazarewicz, T.R. Werner, J.F. 

Berger, C.R. Chinn, and J. Decharge, Phys. Rev. C 53, 

2809 (1996). 

[38] P. MoUer, A.J. Sierk, T. Ichikawa, A. Iwamoto, R. 

Bengtsson, H. Uhrenholt, and S. Aberg, Phys. Rev. C 

79, 064304 (2009). 
[39] R. Smolanczuk, J. Skalski, and A. Sobiczewski, Phys. 

Rev C 52, 1871 (1985). 
[40] Z. Patyk, J. Skalski, A. Sobiczewski and C. Cwiok, Nucl. 

Phys. A 502, 591c (1989). 
[41] A. Dobrowolski, K. Pomorski, and J. Bartel, Phys. Rev. 

C 75, 024613 (2007). 
[42] D.N. Poenaru, M. Miroa, W. Groiner, I. Cata and D. 

Mazilu, Mod. Phys. Lett. A 5, 2101 (1990). 
[43] M. Mirea, L. Tassan-Got, C. Stephan, C.O. Bacri, P. 



Stoica, and R.C. Bobulescu, ,J. Phys. G 31, 1165 (2005). 
[44] M. Mirca, L. Tassan-Got, C. Stephan, C.O. Bacri, and 

R.C. Bobulescu, Europhys. Lett. 73, 705 (2006). 
[45] M. Mirea, L. Tassan-Got, C. Stephan, C.O. Bacri, and 

R.C. Bobulescu, Phys. Rev. C 76, 064608 (2007). 
[46] M. Mirea and L. Tassan-Got, Centr. Eur. J. Phys. 9, 116 

(2011). 

[47] T. Ledorgorber, Z. Paltiel, Z. Fraenkel, and H.C. PauU, 

Nucl. Phys. A 275, 280 (1977). 
[48] A.C. Wahl, At. Dat. Nucl. Dat. Tab. 39, 1 (1988). 
[49] W.J. Swiatecki, and S. Bjornholm, Phys. Rep. 4, 325 

(1972). 

[50] S. Pilat, K. Pomorsky, and A. Staszczak, Z. Phys. A 332, 
259 (1989). 

[51] T. Lodergerbor, Z. Patiel, H.C. Pauli, G. Schutte, Y. 
Yariv, and Z. Fraenkel, Phys. Lett. B 56, 417 (1975). 



